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We calculate in this paper the output power obtainable, versus frequency, 
for TRAPA TT diodes. This high efficiency mode of operation is analyzed 
by means of a simplified model for both germanium and silicon avalanche 
diodes. The model evolved from a study of detailed computer simulations 
of experimental diode-circuit systems. The simplified analysis assumes: 
(i) The avalanche zone transit, phis the recovery time to the swept-out state, 
occurs in a half period of the TRAPATT frequency, (ii) The ratio of 
IMP ATT frequency to TRAPATT frequency is 3:1. {Hi) The diode 
area is chosen to provide 10 ohms negative resistance, a reasonable value 
for microwave circuits. The calculated electrical characteristics agree well 
with experimental observations. Consideration of circuit and thermal 
limitations results in a design for maximum power output for a millimeter 
wave silicon oscillator. Power output in excess of 1 watt CW, with an 
efficiency of 40 percent, is predicted at a frequency of 50 GHz. 

I. INTRODUCTION 

High efficiency oscillations in silicon avalanche diodes were first 
reported in April 1967 by Prager, Chang, and Weisbrod. 1 The new 
high efficiency mode was shown experimentally to be capable of con- 
verting dc to microwave power with efficiency greatly exceeding that 
predicted for the classical IMPATT mode. 2-4 For example pulsed 
bias operation with 60 percent efficiency at about 750 MHz was re- 
ported. 5 Normal IMPATT operation for this structure would be at 
about 4 GHz with perhaps 10 percent efficiency. 
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The first theoretical calculation 8 of the high efficiency mode was 
obtained by computer simulation of carefully characterized experimental 
results 7 on germanium structures. Pulsed bias operation, with output 
power of 7.5 watts at 2.6 GHz with efficiency of 40 percent was re- 
ported. 8 The mode of oscillation was named TRAP ATT, for TRA pped 
Plasma Avalanche Triggered Transit, which describes the new physical 

mechanisms discovered by the computer simulations. 

Characterization of experimental circuits 9 which support TRAP ATT 
oscillations and additional more detailed simulations 1 of the diode- 
circuit system further verify the avalanche-triggered-trapped plasma 
description of the high efficiency mode. 

In this paper a simplified model, previously described, 11 ' 12 for the 
physical processes occurring in TRAP ATT operation is employed in 
order to analyze a wide range of structures and operating conditions. 
Detailed terminal voltage waveforms are calculated, using appropriate 
ionization rates, velocities, and other materials parameters, for germa- 
nium and silicon structures. The avalanche-zone-transit model, first 
described in Ref. 11, is presented in detail, and in a form suited for the 
calculation of the detailed spacial variation of carriers and electric 
field. Analysis of the diode recovery from the "trapped-plasma" state 
follows the treatment presented in Ref. 12. Solutions are obtained by 
applying a terminal current waveform with sharp wavefronts appro- 
priately chosen to represent the effects of the multiresonant circuit to 
the diode terminals. 9,10 A power generating cycle commences with the 
application of a large terminal current to a completely swept-out diode 
structure in which the only carriers available for impact ionization result 
from thermal generation via defect centers. A rapidly moving avalanche 
zone results which sweeps through the depletion region from the p-n 
junction to the heavily doped substrate-base interface. A dense electron- 
hole plasma is left behind the zone and constitutes a high conduction 
current essentially zero voltage state which lasts for approximately a 
half period. Recovery to the essentially zero conduction current-high 
voltage state is also controlled by terminal current waveform. The 
avalanche zone transit occurs in a time considerably less than the transit 
time of a carrier moving at saturated velocity. The sweep-out recovery 
period cannot be less than the transit time. 

The TRAP ATT phenomena described above cannot exist as a small 
signal effect. Instead it results as the natural (and ultimate) large signal 
limit of the chain of events which occur when a sufficiently large drive 
is applied to a swept-out avalanche diode. Obviously, a wide range of 
characteristics result over the range from small signal to large signal, 
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and between single resonance and multiresonance circuits. In this paper, 
we calculate only results in the range of the TRAP ATT large-signal 
limit. Furthermore, to facilitate comparisons between germanium and 
silicon semiconductors, structure designs, bias levels, and operating 
frequencies, we present detailed results for only a particular terminal 
current waveform. This waveform, however, would be expected to 
permit the highest power output at a given frequency, and analysis 
of the output spectrum agrees with experimental circuit characteristics. 
This waveform is a perfect square wave of terminal current chosen so 
that the avalanche zone transient and the subsequent sweep-out of the 
generated electron-hole plasma occurs precisely over the positive half 
period. The negative half period is one of high reverse voltage and low 
saturation current. 

The next two sections (II and III) are mathematical and develop the 
expressions* which are evaluated for a wide range of parameters in 
Section IV. Section II develops in detail the mathematical character- 
ristics of the avalanche zone transient, which occurs when the step in 
current is applied to a completely swept-out diode. Section III is a 
development of the mathematical expressions which describe the re- 
covery period back to the completely swept-out state. Section IV con- 
tains: (i) A summary of the expressions developed in Sections II and 
III, required for the complete analysis. (Readers not interested in 
mathematical detail may omit Sections II and III.) (ii) A discussion of 
the structure design which results as a requirement imposed by the 
particular current waveform. (Hi) Results of efficiency, input power 
density, and diode impedance versus operating frequency for two typical 
structure designs, (iv) For particular optimum structures, results of 
depletion layer thickness, epitaxial impurity concentration, diode junc- 
tion area, average voltage, average current, and output power versus 
operating frequency. Section V presents a discussion of these results, 
and their comparison with experimental information. 

II. SIMPLIFIED MODEL FOR ZONE TRANSIENT 

Selected frames of a computer simulation, 13 shown in Figs. 1-3 
illustrate the zone transit of TRAP ATT operation. Shown are solutions 
for hole and electron concentrations and electric field profiles at suc- 
cessive instants of time. A zone transit begins, as shown in Fig. 1, when 
a large current charges a previously swept-out diode structure until 
impact avalanche generation of carriers begins. For circuit conditions, 

* The results of the analysis are identical to those presented in Ref. 12, although 
the derivation follows from a different point of view. 
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Fig. 1 — Start of zone transient. Hole concentration ( + ), electron concentration 
' — ), electric field (solid) versus distance. 



which maintain this large current, the avalanche zone sweeps through 
the depletion region, as shown in Fig. 2. A dense electron-hole plasma 
is left, as shown in Fig. 3, after the completion of the zone transit. The 
computer simulations utilize a model for the semiconductor diode, 5 
which is quite general and has been quite successful in simulating large- 
signal operation of the various modes of power generation in avalanche 
diodes. Simulated solutions are obtained for the evolution in time of 
the diode variables (hole and electron concentrations and electric fields) 
and circuit variables (terminal current and voltage). The main disad- 
vantage of the program is the cost, which generally prohibits a large- 
scale parameter investigation of the diode-circuit system. 

In this section the expressions are derived which form the basis for a 
large-scale parameter study for the zone transient. The treatment is 
necessarily simplified, and it is assumed that holes and electrons have 
identical properties, except for sign, and in particular the ionization 
rate a and carrier velocity v are identical functions of the electric 
field S. The one-dimensional continuity, transport, and space-charge 
equations are simplified to 
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Fig. 2 — Zone transient period. Indicated are multiplied camera p m and r? m , plasma 
concentrations p p and n p , and plasma field 8 P , resulting from initial concentration 
j t and zone velocity v,. 
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Fig. 3 — Termination of zone transient. 
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t* -#/+«& + «-£. (i) 

fg-f, + «0. + U+fj». (2) 

jv = gpv, (3) 

in = qvn, (4) 

|£ = i (p _ n + Do p) i ( 5 ) 

ox e 

in which diffusion transport has been neglected, g T represents thermal 
generation via defects where significant, and Dop is the net ionized 
impurity concentration. 

A reduction in the number of required equations can be obtained by 
differentiating the space-charge equation and substituting from the 
expression which results when the two continuity equations are added: 

• = (■$- "bP^'**). « 



or 



(„||) = 9 £ (l ,. Dop) + 2 , r + 2a ,_A(i) 



where ; c is the conduction current (j e = ; p + j„) and 

where «/r is terminal current density. 

A single equation in S results by substituting equation (9) into (8), 
but is of no particular advantage. The solution for S and j e is instead 
obtained from the following two equations, where doping gradients 
have been neglected: 

f t = 2v(g T + aj c ) +*£(/*- /.) 



2 a 2 s 

- &) -7-3 

dx 



+ ^g(fDop-.g). (10) 



c ff m **-'•• Cll) 
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Equation (10) is obtained from equation (8) by expressing dv/dx in 
terms of an incremental mobility v' = dv/dS, that is, dv/dx = v' (d&/dx), 
and (dv/dt) = v'(d&/dt). 

2.1 Coordinate Transformation 

The zone transient period of a TRAPATT oscillation is distinctively 
characterized by the rapid movement of an avalanche zone from the 
p-n junction to the base-substrate interface, 11,12 in which computer 
simulations 6 ' 10,13 indicate that all quantities at a time t + A£ have 
profiles related to those existing at time t Q , by the relationship 

S(x, t 4- A*) = &(x - v, At, t a ), 

where the electric field is used for illustration. Incorporation of this 
"experimentally" observed phenomena into the mathematical model 
results in the elimination of all partial derivatives with respect to time, 
and replacing them with an ordinary special derivative. That is by 
making a change of variable z = x — v.- M, 

T f ->-v.£ (12) 

dt az 

where v, is the velocity of the zone and is assumed constant. Incorpora- 
tion of transformation (12) into equations (10) and (11) yields 

-^ = ^.^-2v.(9r + aj c ) (13) 

v, az az 



and 



— & '~7~ = Jt - jc , (14) 



for the case where the velocity is saturated at value v, . The velocity of 
the zone is readily obtained 11,12 from equation (14), when evaluated in 
the region in front of the avalanche, where j, <<C Jt and d&/dz = — 
(q/e)N A : v, = J T /qN A . If thermal generation is neglected, once the 
avalanche zone is formed, equations (13) and (14) can be separated 
and integrated to yield the detailed structure of the avalanche zone 
for any a (8) function, and most interestingly predicts the magnitude 
of multiplication of conduction current through the zone — independent 
of the exact form of the a(&) function: Substitution into equation (13) 
of the derivative of equation (14) with respect to z yields 
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Further substitution of dz from equation (14) yields 



v, v,\ dj r (j r - J 



~ a3r V. vj dZ V ev t 
which can be rearranged as 

9.AQ rfE = (& - ^)(^ I L: ) dj< . (16) 

With the condition that the conduction current in front of the zone 
has constant value j, (as a result of thermal generation during the 
previous clear-out sequence), the relationship between value of electric 
field in the zone and conduction current is given by 



f"«»<-fc-5)l>' h ffi-* + ' 



(2*0- (17) 



From equation (15) we see that v z must be greater than v, in order to 
multiply up conduction current from the value j, . The total current 
multiplication j cm across the entire zone is obtained by evaluating the 
upper limit of the left side of equation (17) also at zero field, in which 
case the left side equals zero and j cm is given by 

3„ m = he <i "- i - /jT \ (IS) 

which shows that the effect of an increase in j, , is a decrease in j cm ." 

2.2 Individual Hole and Electron Concentrations 

The individual hole and electron concentrations can be obtained from 
the value of j c . The expressions are derived for two regions. One is 
anywhere within the zone where v = v, , and the other is to the left of 
the zone where v may be less than v, . 

For v = v„ : 

jc = qvJp + n), 

which can be solved simultaneously with 

obtained by equating equations (5) and (14), as 
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2q \v, v t 

Therefore if the detailed variation of j c is determined, for a particular 
a(£) function, from equation (17), equations (19) and (20) yield the 
detailed profiles of carriers through the zone. In particular the values 
of holes and electrons, see Fig. 2, when j„ = j em are simply 



*-fc 


i L .+i 


- - 1 


\i'„ v.. 



(21) 
(22) 

-(/ V'j v./ 

2.3 For v ^ v, : 

In this region impact ionization can be neglected and we return to 
the individual hole and electron continuity equations (1) and (2) 
written as 

--S--5W. < 23 > 

dz dz 

where v is the carrier drift velocity. Inspection of the above equations 
reveals that any solution must satisfy the following constraints: 

r„ = p(r ; -v), (25) 

C = n(v t +»), (26) 

where c p and c„ are constants. The validity of equations (25) and (26) 
can be checked by differentiation. 

The value of c p and c„ can be obtained by substitution of p m and ??,„ 
from equations (21) and (22), and reveals that c„ = c„ = c : 



The values of p and ?* for any v 2s y, are given by 



(27) 



" = ^-y/ ( "'-"< 6 »' (2S) 



•-*£-?.)/*-•-«»• (29) 
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2.4 Trapped Plasma Solution 

The profiles of holes, electrons and field in the transition region over 
which the carriers change from values p m and n m to their values in the 
trapped plasma, p p and n p , can be obtained from the simultaneous 
solution of equations (28) and (29) and equation (14) for any arbitrary 
v(&) relation. However the solution in the trapped plasma is inde- 
pendent of the exact form of y(8), and the plasma concentration of holes 
and electrons is even independent of the value of low field mobility. 

In the trapped plasma the values of holes, electrons, and field (p p , 
n p , and S p ) are related by 



p v - n p = N A , 

Jt = gv p (p p + n v ), 
v P = n& P . 
The solution for v p is obtained from equations (28), (29) and (30): 



Vr, = 



1c. 



2qN a Xv, v, 



1 + 



2gN n v ! 



• (v, v,\ 



- 1 



(30) 
(31) 
(32) 

(33) 



Once v v is found, the values of p p and n p are given by 



»-*£+' 



n„ = 



N, 



^-1 



(34) 



(35) 



2.5 Discussion of Zone Transient 

A simplified model for the zone transient which results from the 
transformation, equation (12), has been presented. In the numerical 
analysis the calculation begins at the point where the electric field is a 
maximum and j e = J T • Solutions are then obtained by the simul- 
taneous solution of equations (16) and (14) for j e and £ as functions of 
position to the left and to the right of this point. The hole and electron 
profiles are then obtained from equations (19) and (20). The value of 
plasma velocity v p is obtained from a solution of equation (33), and 
the plasma density from equations (34) and (35). 

It is of great significance to note that the solution in the trapped 
plasma is independent of the exact form of the a(S) and «(8) relations 
and is determined uniquely by the value of ;'. , N A , and J T . 

A simplification for the trapped plasma solution results for the condi- 
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tion v p « v, . Equations (33), (34) and (35): 

(36) 



JE — 



»• . v z v. 
X 



n„ = 



fefe-:-:)- 1 )- < 38 > 



In the next section we present an analysis of the recovery period of the 
diode from the trapped plasma state to the state in which the diode is 
completely free (except for thermal generation) of mobile carriers. 

III. RECOVERY TRANSIENT 

In this section the transient solution is calculated for a plasma filled 
n + pp + diode after the application of constant reverse bias current J R , 
which is a good approximation to the recovery portion of a TRAPATT 
oscillation. The analysis also applies to the devices treated by Benda 
and Spenke for the reverse recovery in power rectifiers, 14 and by Moll 
and Hamilton for step recovery diodes. 15 However, the recovery from 
the large density plasma to a swept out state of large voltage, which 
occurs in TRAPATT operation, differs substantially from the earlier 
treatment 14,16 which applies to low density plasmas and slow removal 
rates. 

It is found that the removal rate of the electron-hole plasma and 
build-up of the electric field profile is quite accurately described by a 
simple model 12 which includes the effects of displacement current and 
carrier drift at constant mobility and saturated velocity. For a given 
p-region width the total recovery time can be as small as the transit 
time of the slower carrier if the plasma density is less than a critical 
value and the reverse current equals a critical value. For current 
greater than a critical value avalanche occurs before complete carrier 
removal. The first part of this section pertains to a more general system 
in which holes and electrons are not required to have equal character- 
istics; while the second part is for the special case of equal properties 
and couples the zone transient analysis with the recovery analysis. 

3.1 Part 1 

Figures 4a and 4b are frames from a computer simulation run to test 
the analysis to follow and will serve for notation. 
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Fig. 4 — Recovery transient. Hole concentration, election concentration, and 
electric field versus distance, (a) Plasma velocity period and (b) saturated velocity 
period. 

At t — the diode has a uniform plasma p„ and n p and the center 
region is characterized by: 

Vp = Na + n v i 

Jr = ?£,(m p Pp + MnWp), 

where &„ is the field in the plasma /jl p and n n the hole and elec- 
tron mobilities, and J R the reverse or recovery current. 

The solution for carrier density p, , n, and velocities v L and v% is 
obtained by an integration of the continuity equations over the transi- 
tion region from the plasma to the regions where the field has built up 
to that required for saturated velocity: 



q h! vdx = ~f dj * 



v l and v R are the hole and electron velocities, respectively, in the plasma. 
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and 
or 

(Pv - P,)Vn = VshV, ~ M P p P e„ , 
(>?„ — n a )v L = v SE n, — n„p p & v , 
n p v R = fi„S p n p , 
?y ',. = m p S„p p . 
From the relations above: 

v R = m„S p , (.39) 

r L = M P S P , (40) 

(l'R + Vl) ,.,s 

?,3 = ^(^T^)' (41) 

where Vgn and v SE are the saturated hole and electron velocities. The 
solution as shown in Fig. 4a, with values for v L , v R , p, , n, as given 
above (and verified by computer simulations), is valid until a time t , 
when the plane moving at velocity v L reaches the plane moving at 
velocity v R : 

to - W/(v B + v L ). (43) 

For / ^ /.i the terminal voltage is given by 

V = 2 ^ [vlp s + v\n s + NJ& - v\)]. (44) 

e 1 

For the remainder of the clear-out cycle refer to Fig. 4b. For times 
greater than l u three regions exist: 

Region 1 extends from x = to the plane moving at velocity v SB . 
The field here has slope given by 

d ^= - 5 (A^+0, (45) 

OX e 

and the field & L has a magnitude given by 

«S L /« - J R - qv SE n, = qv L (n, + N A ). (40) 
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Region 3 extending from the plane moving at velocity v s/i to x = W. 
The field here has slope given by 

g-?fr.-*0. (47) 

and the field £« has a magnitude given by 

e& R /t = Jr - qv SIl p, = gv a (p t - N A ). (48) 

Region 2, the remainder, has electric field with slope 

The terminal voltage V for £ ^ £ is obtained, with reference to Fig. 
4b as follows. 

V = (& L + S £ )/2-TF* , 
+ (6* + S„)/2-TF„, 
+ (& E + £//)/2-(TF - TF B - W„), 
or rearranged as 

V = \[Z L W E + S«W'„ + Z E {W - W B ) + &„(W ~ W E )}, (50) 

where 

W E = v L t — v SE (t — t ), 

W H = v R t - v SH (t — t ) 
E E and F,„ are given by equations (46) and (48) and 

£b = Sl - J (#„ + n,)W* . 

£// - S R - 2 (p. - ^)TT ff . 

Equation (50) is valid until W B or W ff equals zero, the time required 
to completely remove electrons or holes. The times required to com- 
pletely clear-out electrons and holes are 

Afr-r"?. (51) 

Vsb "i 

where p, and n, are given by equations (41) and (42). 
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3.2 Minimum Recovery Time 

If the recovery current J R is chosen so that v R = v S b and v L = v SH , 
that is, 

Jr = q(v S nPr, + v SE n v ) = J e , (53) 

the clear-out will proceed at its most rapid rate. Any further increase 
in J R will not speed up the recovery but will result in a dE/dt in the 
plasma which will result in an unnecessarily large value of 8 L or S« near 
the end of the recovery period and premature avalanche before clear- 
out. For J R = J c : 

V. = Vv . 

n, = 7i p . 

A constraint on p p is obtained by requiring the field & L (obtained from 
equation (46) with J R = J e ) to be less than the value obtained there 
in steady-state breakdown S , that is, 

Vv < e&o/gW, (54) 

for t = M H . Since p„ > N A , equation (54) indicates that the diode 
must have a punch-through voltage less than the breakdown voltage. 
The computer simulation shown in Fig. 4 is for silicon with equal 
properties assumed for holes and electrons at 300°K, W = 10 microns, 
p p = 1.52 X 10 16 , n p = 0.514 X 10 16 , N A = 1.0 X 10", and J R = 
1.5 X 10 s A/cm 3 . 

3.3 Part II 

For the remainder of this section it is assumed that holes and 
electrons have equal properties and therefore only one plasma velocity 
v v replaces v R and v L , and one saturated velocity v s replaces v SB and 
v S h in the previous equations. Furthermore, the plasma velocity ob- 
tained during the zone-transit period while current J T flows is not the 
same as the plasma velocity characteristic of the recovery period while 
current J R flows; unless J T = J R . Therefore the following notation 
is used: 

v PT = velocity of plasma while current J T flows, 

v PR = velocity of plasma while current J R flows. 

A factor F is defined by reference to Fig. 5. The breakdown voltage of 
the diode is V , the field at the n + p junction at breakdown is fi , and 
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Fig. 5 — Definition of punch-through factor F. Electric field versus distance at 
breakdown. 



the width of the depletion region is smaller than the value eS, /qN A by 
a factor F > 1. That is F is a measure of the degree of punch-through. 
(In the numerical calculation S and W are obtained by an iterative 
process for any given value of N A and F.) 

It is the purpose of this section to derive design requirements for 
TRAP ATT operation by coupling the zone-transit analysis and the 
recovery analysis. It will be shown that the value of recovery current 
J R which will result in complete removal of the plasma just as the 
diode voltage (and field profile) are the breakdown values V and S is, 

Jh = qN A v s (2F - 1); 

and for the condition that the plasma be extracted in one transit time 
the following constraint must be imposed on the zone transit: 

v t /v PT = 2F - 1. 

The above constraint yields a unique value for J T , for any given 
structure. This value is uniquely determined by material (Ge or Si), 
N A , F, and J s • For an additional constraint, for example that we 
minimize the sum of zone-transit time and recovery time or that this 
time period have a fixed relationship to the IMPATT transit time, a 
unique value of F can be determined. By such constraints then a unique 
(except for weak dependence on J s ) power- versus-frequency curve can 
be obtained for a given material. The derivation of the expressions 
required follows. 
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At the end of the zone transit the diode has plasma density given by 
equations (34) and (35). The clear-out densities p, and n, are 

'■ - N i~ + l ) i^+7 • (55) 

n 8 = A^(^-l)-^~- (56) 

Utilizing equations (55) and (52) the electric field S L , obtained from 
equation (4G), just as the carriers are cleared out has value 



f,,(A/) - qv PR (n. + tf*)(^&) 



and when set equal to S„ (from Fig. 5 £„ = (q/e)N A FW) yields 

(n. + N A ) V -^ = N A F. 

Further substitution from equations (41) and (56) yields 

2F = i + & fca 

or utilizing equations (41) and (42) 

/« - giV^,(2F - 1). (57) 

The equation above results solely from the requirement that right 
when the last remnants of the plasma are leaving the depletion region, 
the diode is at breakdown voltage V . We have the interesting result 
that the degree of punch-through (F) uniquely determines— an upper 
limit on the reverse recovery current, independent of the zone transit 
current J r , ionization rate, mobility or saturation current. 

The critical current for removal of the plasma in one transit time is 

J c = q?.(pv + «p)- 
For p„ and u„ given by equations (34) and (35), 

J c = qv,N A v./v PT , 
and for the condition J R = J c , 

■^- = 2F-1. (58) 

VpT 

Equation (58) imposes a constraint on the zone transit process so that 
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the resulting plasma can be removed in one transit time. From the 
zone- transit analysis, in particular equation (33), a relationship was 
derived relating v t , v, , j cm , N A and v PT . Therefore, combining equa- 
tions (33) and (58) yields a relationship between F, N A , v, , J T and j, , 
which must be satisfied if the plasma is to be extracted in one transit 
time. However, the condition for maximum frequency is one which 
minimizes the sum of three characteristic times: the charging time t x , 
the zone transit time t 2 and the recovery time (compare Fig. 6). 

The power output of an avalanche diode is limited by realizable 
circuit impedance or thermal limitations. Both have the effect of 
limiting the area of the diode. For a given upper limit on circuit im- 
pedance the best design is one which maximizes the product of power 
output times frequency squared. This figure of merit (F.M.) will have 
a maximum value determined by a judicious choice for J T , J R , the 
duration of time they are applied, and the duration of time the break- 
down voltage-saturation current state persists. Such optimum condi- 
tions are not sought in the present analysis, inasmuch as after they 
were found it would probably be quite difficult to design a realistic 
circuit which could provide such a waveform. Instead, optimum values 
of F and current amplitude are found such that the F.M. is a maxi- 
mum for perfect square wave current drive. (This is experimentally 
observed 10 where the square wave is made up of a limited number of 
odd harmonics; the circuit is an open circuit at all even harmonics, 




0.3 0.4 0.5 

TIME IN NANOSECONDS 



Fig. 6 — Example of calculated voltage (solid) and assumed diode current (dashed) 
versus time. Zone transit (v r ), plasma velocity recovery (v p ), and saturated velocity 
(v.), periods indicated. 
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the fundamental has impedance as calculated herein, and the third 
harmonic corresponds to the IMPATT frequency.) The next section 
reports results of a calculation for such a power-frequency characteristic. 

IV. POWER OUTPUT AND STRUCTURE DESIGN 

The application of terminal current J T to a swept-out diode results 
in the zone transients shown in Figs. 1-3. In Section II we presented 
a considerably simplified analysis for the zone transient, as shown for 
example in Fig. 2. In Section III we analyzed the recovery transient, 
as shown in Figs. 4. A recovery current J R = qN A v.(2F — 1), equa- 
tion (57), is required to completely remove the plasma in as short a 
time as possible, such that the diode voltage and field profile have 
just attained the breakdown values. The structure factor F and break- 
down values are defined by Fig. 5. In the first part of this section 
results are presented on a per unit area basis. The voltage waveform, 
as shown for example in Fig. 6, is calculated for various values of N A , 
ranging from 10 15 cm -3 to 10 17 cm -3 , and values of F equal to 1.5 and 
2.0. The current density J T is determined by equation (57). The re- 
sulting waveform is Fourier analyzed, the frequency determined, and 
results presented in Figs. 7-9. In the second part, results are presented 
for an optimally chosen structure: (i) The value of F is found for each 
value of N A such that the ratio of IMPATT frequency (tt transit 
angle) to TRAP ATT frequency is 3 : 1. (it) The diode area is chosen 
to provide 10 ohms negative resistance, a reasonable value for micro- 
wave circuits. These results are presented in Figs. 10-15. 

4.1 Parti 

As shown in Section II the zone transient solution is logarithmically 
related to the thermally generated current j, . In the calculations the 
current j, is given by 

j. = griiW/T. (59) 

The intrinsic carrier concentration is evaluated at 300°C. for silicon 
and 100°C. for germanium, and the lifetime t is 100 nsec. Effects of 
variations of j. on device design are considered elsewhere. 12 The satu- 
rated velocity v, is 8.9 X 10 8 cm/sec for silicon and 6.5 X 10° cm/sec 
for germanium. The ionization rate a = ae~ b/E has a and 6 1.6 X 10 
cm" 1 , 1.9 X 10° V/cm for silicon and 3.5 X 10 5 cm -1 , 7.5 X 10 6 V/cm 
for germanium. 16 

The waveform of diode voltage for a particular structure is shown 
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Fig. 7 — Efficiency versus frequency, F = 1.5 and 2.0. 
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Fig. S — Input (average) power density versus frequency, F = 1.5 and 2.0. 
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Fig. 9 — Diode impedance versus frequency, F = 1.5. 
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Mg. 1(» — Depletion layer thickness \V versus frequency. 
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Fig. 11 — Impurity concentration N A versus frequency. 
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Fig. 12 — Diode area versuB frequency. 
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Fig. 13 — Average voltage versus frequency. 
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Fig. 15 — Output power versus frequency. 



in Fig. 6. The time t x is a depletion layer capacitance charging time, 
required to overdrive the diode to the point at which the avalanche 
zone transient begins. The zone transient lasts until time t 2 . The re- 
covery transient characterized by the plasma velocity v p lasts until 
time t 3 . The remainder of the recovery time, to t 4 , is characterized by 
carriers moving at saturated velocity. After complete sweep-out of 
carriers, the diode voltage equals the breakdown value and the diode 
current equals the saturation value j. , for a half period of duration t t . 
The voltage waveforms are calculated according to the simplified 
models presented in Sections II and III. The voltage for times between 
U and t 2 , the zone transient, are calculated by "sliding" the field pro- 
file, shown for example in Fig. 2, from the point where the maximum 
field occurs at the N + P junction, corresponding to time t t , to the right 
at velocity v, . The time t 2 corresponds to the time when the location 
of the right-most point of the plasma region would enter the substrate. 
The voltage is the area under the field curve, within the confines of the 
depletion region. The voltage thus calculated is only an approximation 
at the beginning of this transient. However, the value of trapped 
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plasma left in the depletion region at time t 2 is independent of the actual 
avalanching conditions, and this plasma constitutes the initial condi- 
tions for the recovery voltage waveforms from time t 2 until t 4 . The 
voltage waveform in the neighborhood of time t x is determined by 
complicated space charge effects shown in Fig. 1. It is felt that detail 
occurring over such a short period of time would not significantly 
affect the overall power output and efficiency calculation, but this is 
one of the major limitations of this simplified analysis. 

Results for two values of punch-through factor F are shown in Figs. 
7 and 8. The frequency is calculated (for particular values of N A ) and 
results shown versus frequency. Figure 7 is a plot of calculated power 
conversion efficiency versus frequency; Fig. 8, dc power density versus 
output frequency. The diode impedance, as calculated by Fourier 
analysis of the waveforms, is shown in Fig. 9 for F = 1.5. 

4.2 Part II 

Results shown in Figs. 7-9 are on a per unit area basis, and no opti- 
mization for the F factor has been attempted. The area of an avalanche 
oscillator is determined by circuit impedance limitations and thermal 
dissipation limitations for CW operation. Results of the calculation 
are presented in Figs. 10-15 for an area chosen such that 10 ohms 
negative resistance at the TRAPATT frequency is obtained. In addi- 
tion, the F factor was chosen by an iterative process for optimization 
of the diode structure. Initially optimization was sought for the cri- 
terion of maximization of the products of output power times operating 
frequency squared. Results clearly indicated that the maximum was 
for values of F only slightly larger than 1.0. For such cases the recovery 
time is very nearly one transit time, the minimum value. The total 
half period of the oscillation would be only slightly larger than two 
transit times, and therefore a ratio of IMPATT-to-TRAPATT fre- 
quency of 2 : 1 could be obtained. Instead results for a frequency ratio 
of 3 : 1 are presented. This will be somewhat less than the optimum, 
but more realistic, as this mode of operation is more easily obtained 
experimentally. 

Device depletion layer thickness W, epitaxial impurity concentra- 
tion N A , and device area A, calculated for the 3 : 1 mode, are shown 
in Figs. 10-12, versus output frequency. Operating characteristics are 
presented in Figs. 13-15. Plotted versus frequency are average voltage, 
average current, and output power. As discussed previously, device 
area was chosen to provide 10 ohms negative resistance at the operating 
(TRAPATT) frequency. No other circuit limitations on area are in- 
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eluded, such as might be imposed by IMPATT considerations. Such 
considerations will be discussed later. The output power calculated for 
a silicon device, if the area is chosen for CW operation, is shown as a 
dashed line in Fig. 15. The choice of maximum device area permitting 
CW operation was obtained from silicon IMPATT diode thermal 
resistance and burnout analysis collected over the past several years 
at this laboratory. 17 Based on this analysis, circuit impedance limita- 
tions determine device area for operating frequencies above about 10 
GHz. In addition, device area may have to be reduced to provide 
sufficient negative resistance to overcome losses at the IMPATT fre- 
quency. This effect has not been considered in the present analysis. 

v. DISCUSSION 

We have presented calculated operating characteristics of the 
TRAP ATT mode of high efficiency power generation in germanium 
and silicon avalanche diodes. The calculated results were obtained by 
analysis of solutions of a simplified model for the avalanche zone 
transit and reverse recovery periods which characterize the TRAP ATT 
mode. The calculated solutions agree well with solutions obtained by 
detailed computer simulation studies, except for detail missing at the 
beginning of the zone transit, as the avalanche zone is forming. 

Characteristics were presented, on a per unit area basis, for diodes 
with two degrees of punch-through, of the electric field at breakdown, 
into the terminating substrate. The silicon structures are more efficient 
than the germanium but require considerably more input power. For 
both materials input power requirements increase with F, the degree of 
punch-through, at a particular frequency. 

A particular design was obtained by choosing F such that the 
TRAPATT frequency would be one-third the IMPATT (tt transit 
angle) frequency of the structure. The calculated epitaxial impurity 
concentration and depletion layer width were presented. The calcu- 
lated results are in excellent agreement with our experimental Ge 
structure (/ = 2.4 GHz, W = 4.5 M m, N m = 4.2 X 10 15 cm -3 ). 8 ' 10 Cal- 
culated results of power output versus frequency were presented for 
diode areas chosen to yield 10 ohms negative resistance at the 
TRAPATT frequency, and for silicon diodes with areas determined 
by thermal limitations. Power output in excess of 1 watt CW, with 
efficiency of 40 percent is predicted at a frequency of 50 GHz. 
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